Dynamic scaling in stick-slip friction 



in 
o 
o 

(N 

S-H 

Oh- 
< 



o 

< 



> 

m 
o 
^t- 
o 
in 
o 



> 

x 



Jens Feder0 Hakon Nordhagen, and Wesley Andres WattersQ 
Physics of Geological Processes, University of Oslo, Box 1048 Blindern, N-0316 Oslo, Norway 

(Dated: February 8, 2008) 

We introduce a generalized homogeneous function to describe the joint probability density for 
magnitude and duration of events in self-organized critical systems (SOC). It follows that the cumu- 
lative distributions of magnitude and of duration are power- laws with exponents a and r respectively. 
A power-law relates duration and magnitude (exponent 7) on the average. The exponents satisfy 
the dynamic scaling relation a — jr. The exponents classify SOC systems into universality classes 
that do not depend on microscopic details provided that both a < 1 and r < 1. We also present new 
experimental results on the stick-slip motion of a sandpaper slowly pulled across a carpet that are 
consistent with our criteria for SOC systems. Our experiments, as well as experiments by others, 
satisfy our dynamic scaling relation. We discuss the relevance of our results to earthquake statistics. 

PACS numbers: 05.65.+b,89.75.Da,46.55.+,91.30.-f 



INTRODUCTION 

The force required to pull a block along a surface is pro- 
portional to its weight and independent of contact area. 
The sliding motion of the block depends on the elastic 
properties of the string pulling the block, the block itself, 
and the surface, and takes many forms Q- S Q| • F° r a 
stiff system, periodic slips with a velocity dependent am- 
plitude have been observed 0, . Intermittent stick-slip 
friction has been observed in experiments where a highly 
compliant string was used to pull a piece of sandpaper 
across a soft carpet [6j. We extend these measurements 
to resolve the slip duration and find that our system is 
an example of a self-organized critical (SOC) system 0. 

We propose a testable criterion for scale-free SOC sys- 
tems: The joint probability density of slip magnitude 
and duration is a generalized homogeneous function, the 
exponents of the cumulative magnitude and duration dis- 
tributions must be less than unity, so that moments of 
the distributions do not exist, and a dynamic scaling re- 
lation between the exponents must be satisfied. It follows 
that periodic stick-slip motion 0,O| cannot be SOC. The 
experimental results presented here satisfy the definition 
of SOC accurately. 

The scaling relation is also satisfied by the magnitude 
and duration exponents of rain showers [8j, and by ex- 
ponents for slips in a granular medium [9|. In our view, 
earthquakes are not SOC since the currently accepted 
exponents are larger than one. 



EXPERIMENTS 

A weighted sandpaper was pulled by an elastic nylon 
line slowly (9 /xm/s) across a carpet [lOj. The sandpaper 
moved by a stick-slip process. We measured the force re- 
quired to hold the carpet in place (see Fig.^l at intervals 
of St = (lfflz)^ 1 , that is, every 0.6 /im on the average. 
For each slip, we observed an abrupt decrease in force, 



AF, and for the first time we have been able to measure 
the duration, t, of slips. 

The force per unit area (stress), f(r), between the 
sandpaper and the carpet is a function of position, r, 
within the surface area, A, of the sandpaper. The in- 
tegral of the vertical component of f(r) is a constant 
equal to the weight of the sandpaper slider (mass m): 
mg = j A f(r)-n d 2 r, where g is the acceleration of grav- 
ity, and n is the surface normal. In a slip, the contact 
between the sandpaper and the carpet is lost, or reduced, 
over an area S < A, that is not necessarily connected. 
The force changes by 5f(r) — /1 u(r)/h, where u(r) is 
the local slip vector, in the plane of the carpet, h is the 
carpet thickness, fi is the shear modulus of the carpet, 
and u(r)/h is the slip-associated change in shear strain 
at r. The slip leads to an elastic increase in shear strain, 
u e (r)/h outside S, and increased elastic stress takes more 
of the external load. Integrating Sf(r) over the area of 
the sandpaper, A, we find that 

m = hAF = h Sf(r) ■ d 2 r 

J A 

= [i I u{r) • d 2 r + /1 / u e (r) • d 2 r 
Js JA-S 

= v'S(u) s ■ 

Experimentally we measure the first integral. The next 
integrals represent a model where slip occurs over an area 
S, and an increased elastic strain elsewhere in A. The 
integral over A — S is proportional to the integral over S, 
and {u)s is the mean slip averaged over the slip area S. 
The last expression is the standard form of the (scalar) 
seismic moment Mo with an effective shear modulus // 
[ill . We conclude that by measuring A_F we obtain 
the seismic moment, m, of the slip directly. 

A slip starts at t\ if the force difference between two 
consecutive measurements SF(ti) = F(ti + t) — F(ti), 
changes sign from positive to negative at t\, and ends 
at £2 = t\ + t when SFfo) becomes positive. A slip is 
therefore characterized by the decrease in force AF = 
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FIG. 1: Force as a function of time from an experiment. The 
initial rise and drop in force at the beginning of the exper- 
iment is characteristic of the interval of time in which the 
system "forgets" its initial state. Later slips of size AF and 
duration t, are studied here. 
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FIG. 2: The cumulative distribution, P(M > m), as a func- 
tion of the seismic moment, m. The observed cumulative dis- 
tribution is consistent with a power-law, P(M > m) ~ m~ a , 
with a = 0.67 ± 0.03. The straight line has a slope of —0.67. 



F(ti) — F(t<z) : that is, by the seismic moment m = hAF 
and duration t = t% — 1\. We found that this procedure 
gave the same scaling exponents as more complicated 
ways of estimating slip magnitude and duration, such 
as sliding average and other filtering procedures, that 
we have tested. Of course, each slip consists of events 
whereby fibers lose their grip on the sandpaper one by 
one; the sound from slips was easy to hear. From this 
point of view the minimum slip consist of one single fiber 
losing its grip. Wc have seen that there are bursts of such 
micro-slips in the slips identified by us. There was some 
very low activity even in the intervals between slips. Of 
course, electrical noise from the strain-gauge, and mea- 
surement system limits the resolution so that we cannot 
separate the smallest slips from electrical noise. For the 
carpet used in the experiments reported here, in no case 
did all fibers lose their grip on the sandpaper, so that 



FIG. 3: The cumulative distribution, P(T > t), as a function 
of slip duration t. The observed cumulative distribution is 
consistent with a power-law, P(T > t) ~ t~ T , with r = 0.96 ± 
0.03. The straight line has slope —0.96. 



the sandpaper was never displaced more than one fiber 
length with respect to the carpet. In the case of carpets 
for which this did sometimes occur, a power-law distri- 
bution of magnitudes and durations was not observed. 

Figure [21 shows that the cumulative distribution of 
magnitude, averaged over four independent experiments 
with a total of more than 50 000 slips, is consistent with 
the scaling form 



P(M > to) ~ to" 



a = 0.67 ±0.03 



(1) 



The deviations from the power-law fit (dashed line seg- 
ments were not part of the fit) at large m are due to 
finite-size effects, and the increase at small m is due to 
noise, which masquerades as small slips. Similar devia- 
tions for small and large values are seen in Figs. 01 and 01 
Since power-law distributions are stable under addition 
[lij , the effect of noise and finite sampling rate does not 
affect the slope of the distributions, but reduces the range 
over which scaling may be observed. 

Figure |3| shows that the cumulative distribution of 
event duration time t is consistent with the scaling form 



P{T > t) ~ r 



0.96 ±0.03 



(2) 



We also found that slips with a large seismic moment 
lasted longer than small slips. Figure 01 shows the esti- 
mates of most probable values (mo, to) of the joint prob- 
ability density p(m, t) as points, and estimates of the 
standard deviations (for bins that contain more than six 
measurements) around these points as the shaded area, 
that is, we have estimated to to be the average duration 
given m = mo- 

The results shown in Fig. 0] are consistent with the 
power-law form relating the most probable pairs (toq, to): 
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for m > 0.15 dyn cm. 



7 = 0.69 ±0.03 , 
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FIG. 4: The joint probability density, p(m, t), of seismic 
moment, m, and slip duration, t, is illustrated by showing 
the most probable values (mo, to) as points. The shaded area 
represents (m, t) values within one standard deviation around 



(mo, to)- The straight line is a fit of to 
7 = 0.69. 



mjj with a slope 



THEORY 

SOC, as originally denned Q, requires power-law scal- 
ing both in magnitude and duration of the events. We 
show that assuming the joint probability density p(m, t) 
to be a generalized homogeneous function, it follows that 
the magnitude and duration distributions are power-laws, 
and average duration as a function of magnitude is also 
a power-law. The measured exponents satisfy a scaling 
relation that follows from the assumption that the joint 
probability density, p(m, t), is a generalized homogeneous 
function, that is, for any positive value of A: 



p(m, t)dm dt = A a p 



to t 
A~' A 7 



(4) 



Here A is a number that represents the change of scale 
provided to and t are given in dimensionless units. Such 
functions have been discussed extensively in the context 
of scaling at the critical point of second order phase tran- 
sitions |lq |. However, for non-equilibrium systems, of the 
kind discussed here, the theoretical basis for Q has not 
been established and it thus remains a proposal. 

The probability density, p(m), for slip magnitude is ob- 
tained from eq. (@J by integrating over the slip duration 



p(m)dm = / p(m, t)dm dt 



it-o 
= X-^+^dm 



P 



v=0 



(?-) 



dv 



where we have introduced the integration variable v = 
i/A 7 . Since A may be chosen to have any positive value 
we set A = to, and find 



p(m) = m 



= rr ) -( 1 +«) / p(l,v)dv ~m-( 1+a K 
v=0 



Here the integral is a constant, and it follows that the 
cumulative distribution, P(M > to) = J p(M) dM ~ 
to - ", has the scaling form given in eq. Q. 

The probability density, p(t), for the slip duration t is 
found, in a similar way, by integrating eq. over to: 

/>oo 

p (t) = t-( 1+a '^ / P (v, i) dv ~ r (1+Q/7) . 



Again the integral is a constant, and it follows that the 
cumulative distribution, P(T > t) = j t °° p{T) dT ~ 
t~ a h ^ nas the scaling form given in eq. (j2J with r = a/7, 
that is, we have the dynamic scaling relation 



a = 7T 



(5) 



This relation also follows by assuming the scaling forms 
JU ® and requiring consistency for the conditional ex- 
pectation values [l6j| . 

Equation JHJl follows from our proposition Let 
(toq,^) be a (to, t) pair where p{m,t) has a local max- 
imum corresponding to the mostprobable value for m 
given t = £q or t given to = to,q. Then, for the generalized 
homogeneous function J3J|, m = TOq/A and to = ig/A 7 
also corresponds to a local maximum in the joint proba- 
bility density for any positive value of A; it follows that 

1/A 7 - to/*o - ("VmS) 7 , 
which is indeed eq. 10. 



to be dimensionless, we may chose both TOg and tj$ 



Since both to and t are taken 

to be 

unity for the purposes of the scaling relations. Equation 
states that the joint probability density is invariant 
with respect to the affine transformation to — + to/ A, t — > 
t/X 1 , and p — > p/X a ; it follows that any expectation and 
other characteristic, such as the full width at half height, 
will scale and have exponents that only depend on a and 
7. It follows that the width of the distribution in Fig. 01 
should be constant in the log-log plot. 

If there is a characteristic magnitude, or duration, of 
slip events, then any model of the process must include 
the characteristic values. However, if a < 1 and r < 1, 
no characteristic value exists. Consider the expectation 
value of the magnitude 



(to) 



mp(m)dm 



to "dm 
a < 1 and mj, — > 00 , 



1 — a 



a > 1 and m f 







(6) 



Here mj is the maximum event magnitude, which di- 
verges with system size L, and m e is the smallest, "atomic 
scale", event size. If a < 1, then (to) fails to exist and 
diverges with system size (note that (m) — » Iiito^ as 
a — > Similarly, if r < 1 then (t) — > 00 as sys- 

tem size increases and the system is also scale-free in the 



4 



time domain. In this case the system is scale-free and has 
no characteristic event magnitude or duration; and any 
model of the system must also be scale-free and cannot 
depend on atomistic details. We consider such models 
not to be SOC. 

We therefore consider systems to be self-organized crit- 
ical provided that eqs. J3J and (J5J are satisfied, and the 
first and higher moments of both the magnitude- and the 
duration probability densities diverge with system size. 

From the modeling point of view, we expect that sys- 
tems with finite first moments, that is, a > 1, and/or 
r > 1, require "atomistic" details since the distributions 
then are dominated by the small scale and/or short time 
behavior, i.e., dominated by m e in © even if mi — > oo. 



TABLE I: Scaling exponents for magnitude and duration 



Distribution 




Carpet Rain Granular Earthquakes 

® [9] \n,u} 


P(M > m) ~ mr a 


Q 


0.67 


0.4 


0.94 


2/3 


P(T > t) ~ t~ T 


r 


0.96 


0.6 


1.08 


2 




7 


0.69 


0.7 


0.87 


1/3 


a — 7T 


7T 


0.66 


0.4 


0.94 


2/3 



DISCUSSION 



Mq [llj, which has the Gutenberg-Richter distribution 
P(M > M ) ~ Mq B with the exponent B — a = 2/3 
[ill ll2T |. Implicitly it is assumed that Mq is the only 
stochastic variable for earthquakes, other quantities such 
as duration t, and the area S are "mechanically" related 
to Mq. The duration is t = L/V, where V is the rupture 
velocity, and L — S 1 / 2 ; also M oc 5 3 / 2 . These scaling 
relations are consistent with observations ^ follows 

1 /3 

that t oc M ' , i.e., j = 1/3; and r = 2. Consequently, in 
our view, earthquakes are not SOC since they are not de- 
scribed by a joint probability density for magnitude and 
duration that is a generalized homogeneous function, and 
because r > 1. Models of earthquakes must include "mi- 
croscopic" details for the duration, and be scale free for 
magnitudes, in contradiction with the mechanistic rela- 
tion t = M 1/3 . 

The original paper Q on SOC introduced a cellular 
automaton model for avalanches (with a — for a (50) 2 
system). However, the exponents are not well determined 
since they depend strongly on system size |2(ij |. For a 
(2048) 2 system a = 0.13, 7 = 0.61, and r = 0.20 consis- 
tent with eq. JSJ. The conservative OFC model [ilj has 
exponents a = 0.25, 7 = 0.42, and r = 0.52 satisfying 
the scaling relation. Here duration was taken to be the 
number of cellular automaton updates required to end 
an avalanche. Both models are SOC with our definition. 
The models just discussed cannot be in the universality 
class of our stick-slip motion, or earthquakes, since they 
have too small values for a. 



Table [I] collects values for the scaling exponents and 
tests the validity of the dynamic scaling relation eq. JSJ. 
In our experiments we have measured the exponents a, 
t, and 7. We find that both a < 1 and r < 1, and that 
the scaling relation a — 7T holds within experimental 
accuracy. We conclude that the stick-slip process studied 
is SOC. 

The dynamic scaling relation may also be tested on 
other "avalanche phenomena" where size and duration 
are known to be power-law distributed such as the size 
and duration of rain showers || , and the slips in a gran- 
ular medium under shear [jj. We note that the granu- 
lar slip experiments have r = 1.08 and the system 
is therefore not scale free unless r < 1, which is within 
experimental error. 

The exponents take different values for different sys- 
tems. In the language of critical phenomena ^^|, we 
would say that stick-slip motion and rain showers belong 
to different universality classes, since the exponents dif- 
fer. We expect that for each SOC universality class there 
exist a scale free dynamic model characterizing that class. 

Earthquakes have long been recognized as resulting 
from stick-slip friction [171. Earthquakes may repre- 
sent an example of SOC 0, E, El- The magnitude of 
earthquakes is characterized by their seismic moments 



CONCLUSION 

We conclude that stick-slip motion in the system we 
have studied exhibits power-law scaling that is well de- 
scribed by a joint probability density for magnitude and 
duration that has the form of a generalized homogeneous 
function. We derived a scaling relation between the dy- 
namic exponents r and 7, and the magnitude exponent a 
that takes the simple form a = 77". We have found other 
examples in natural and experimental systems that sat- 
isfy this scaling relation. We show that our analysis leads 
to open questions for earthquakes. 

We have proposed that self-organized critical (SOC) 
systems are characterized by 

• a joint probability density for magnitude and du- 
ration that is a generalized homogeneous function, 

• with exponents a < 1 and r < 1, and 

• that satisfy the scaling relation a = 77". 

We thank Kim Christensen and Y. Kagan for discus- 
sions. WAW thanks the Fulbright foundation for a fel- 
lowship. 



5 



* Electronic address: feder@fys.uio.no 

* Also at: Department of Earth Atmospheric and Plane- 
tary Sciences, MIT, Cambridge, MA 02139, U.S.A. 

[1] F. P. Bowden and D. Tabor, The Friction and Lubrica- 
tion of Solids. Part I & II (Oxford University Press, 1954 
& 1968). 

[2] B. N. J. Persson, Sliding Friction: Physical Principles 
and Applications ( Springer- Verlag, Heidelberg, 1998). 

[3] E. Gerde and M. Marder, Nature 413, 285 (2001). 

[4] T. Baumberger, F. Heslot, and B. Perrin, Nature 367, 
544 (1994). 

[5] F. Heslot, T. Baumberger, B. Perrin, B. Caroli, and 

C. Caroli, Phys. Rev. E 49, 4973 (1994). 
[6] H. J. S. Feder and J. Feder, Phys. Rev. Lett. 66, 2669 

(1991), erratum: Phys. Rev. Lett. 67:283, 1991. 
[7] P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 

59, 381 (1987). 
[8] O. Peters and K. Christensen, Phys. Rev. E 66, 036120 

(2002). 

[9] F. Dalton and D. Corcoran, Phys. Rev. E 63, 061312 
(2001). 

[10] The setup is similar to that of |6j, but the pulling speed 
was only 9 /im/s. In a typical experiment, a circular sand- 
paper (60 grit), 12.5 cm in diameter, was mounted on a 
heavy stiff support (total weight 630 g) and pulled by an 
elastic nylon fishing line (1 m long and 0.6 mm in diam- 
eter, with an elastic constant of 420 N/m) A new carpet 
was used for each experiment. The carpet consisted of 
wool and synthetic fiber bundles in 4 mm high loops at- 
tached to a 1 mm cloth layer backed by a 3 mm latex 



layer. Approximately 1300 fiber bundles were covered by 
the sand paper. For each fiber bundle, 10-30 fibers came 
into contact with the sandpaper. That is, the number of 
fibers in contact with the sandpaper was in the range 
15 000-45 000. The force on the carpet was measured us- 
ing a strain gauge (Omega LCL-010 mounted between 
the movable aluminum plate, supporting the carpet, and 
the fixed motor support) at 14 Hz with a Keithley 2002 
digital volt meter. 

[11] C. H. Scholz, The Mechanics of Earthquakes and Faulting 
(Cambridge University Press, 2002). 

[12] H. Kanamori and E. E. Brodsky, Rep. Prog. Phys. 67, 
1429 (2004). 

[13] In reference |g the average speed was 270 /im/s, com- 
pared with only 9/mi/s here. The magnitude exponent 
was q = 0.79 ±0.05, which is above, but consistent with, 
the value found in the extensive experiments presented 
here. 

[14] B. B. Mandelbrot, The Fractal geometry of Nature (W. 

H. Freeman & Co, 1982), chap. 38 and pp. 367-. 
[15] C. Domb, The critical point (Taylor & Francis Ltd, 1996) . 
[16] K. Christensen, H. C. Fogedby, and H. J. Jensen, J. Stat. 

Phys. 63, 653 (1991). 
[17] C. H. Scholz, Nature 391, 37 (1998). 
[18] P. Bak and C. Tang, J. Geophys. Res. 83, 15 635 (1989). 
[19] P. Bak, K. Christensen, L. Danon, and T. Scanlon, Phys. 

Rev. Lett. 88, 178501 (2002). 
[20] K. Christensen and Z. Olami, Phys. Rev. E 48, 3361 

(1993). 

[21] Z. Olami, H. J. S. Feder, and K. Christensen, Phys. Rev. 
Lett. 68, 1244 (1992). 



